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STRESS-STRAIN AND FAILURE PROPERTIES 
OF GRAPHITE/EPOXY LAMINATES 

ABSTRACT 

The results of a series of tensile tests on qraphl te/epoxy [0°]g S 
and [0°/+30°/0°]2^ laminates at rates varying from 0.0C2 in/mln to 2 In/ 
min are reported. The loads are applied at various angles to the fiber 
directions In each case. The rate dependent behavior of the stress- 
strain response Is assessed. Evidence Is presented to indicate that 
failure first occurs on Inner plies. Also, evidence is presented to in- 
dicate that. In some cases, moduli increase with increased stress or 
strain level. Lamination theory is used to predict moduli and compari- 
sons with experiment are given. Also, lamination theory is used in con- 
junction with three failure theories to predict ultimate strenqths with 
varying degrees of success. 
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Introduction 

The proper characterization of advanced laminated composite 

materials Is quite complex due to their Inhomogeneous and anisotropic 

nature and due to the coupling effects caused by the lamination process. 

The latter often results In delamlnatlons or stacking sequence effects 

1 ? 

as shown by Daniel, et al. Numerous testing techniques and analytical 
techniques-* have been proposed for the characterization of composite 
materials. Basically, the behavior of composites can be assessed from a 
micromechanical or a macromechanical viewpoint. Essentially, the rule of 
mixtures Is representative of the former approach and lamination theory 
Is representative of the latter. 

The macroscopic approach of lamination theory has been used In the 
present effort. In brief, the goal of this study was to Investigate the 
tensile properties of [0°]ss and [0 o /±30°/0°]25 graphite/epoxy laminates. 
This Included moduli and strengths and their prediction by appropriately 
applying lamination theory In conjunction with several failure theories. 
Further, the viscoelastic or rate dependent behavior of the laminates 
Investigated (Hercules AS fiber and 3501 resin) was to be determined. 

Analytical Considerations 

The constitutive relations using classical lamination theory for an 
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orthotropic lamina In a state of plane stress can be written as 
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The C^j are the stiffnesses and K refers to the lamina. Only four 

elastic constants are required, namely, , Gl2 an d v 12 w ^ iere 

subscripts 1 and 2 refer to the fiber direction and transverse to the 

fiber direction respectively. The latter constants are related to the 
3 4 

stiffnesses (C^ j ) . * The 1-2 coordinate system Is along the axes of 
material symmetry. This local coordinate system can be transformed Into 
the global laminate x-y coordinate system which gives the followlnq 
constitutive relations for a lamina 



C 11 C 12 C 16 
Ci 2 C 22 C 26 
Cl 6 c 26 C 66 



( 2 ) 


Here the Em are the transformed or reduced stiffnesses. 

The laminate stiffnesses, A^j, can be calculated from the trans- 
formed laminae stiffnesses, fij. For a symmetric laminate, the In- 
plane loads are related to the In-plane laminate strains, {rj}, by 


{Nj } = C A ij] 1 and j * x or y (3) 


where {N^} are the resultant forces per unit length of the laminate. 

(For uniaxial tension N x f 0 and N y = N xy * 0.) Thus, the In-plane 
strains can be calculated knowing the In-plane laminate stress resultants 
and the properties of a single lamina. Further, this information can be 
used to calculate the stresses In any layer or lamina In terms of the 
externally applied loads. Since Interlaminar effects are not included, 
the analysis predicts Identical response for a given laminate regardless 
of stacking sequence of individual plys. 
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Numerous failure criteria are presently available which can be 
classified as either having Independent or dependent failure modes. The 
maximum stress or maximum strain criterion are ones with Independent 
failure modes. The criteria proposed by Ashkenazi,^ HI 1 1 ^ and Puppo and 

O 

Evensen are examples of failure theories with dependent failure modes. 

Between the two classes, It would seem that those with dependent failure 

modes would be more appropriate for comoosite materials. Experimental 

a 

results tend to confirm this assumption. 

Ashkenazi's theory assumes a macroscopic continuum In which the 
strength properties are fourth order tensors and In which environmental 
effects are neglected. His criterion for a plane orthotropic material 
under uniaxial load In an arbitrary direction Is 
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where o x is the applied normal stress making an angle e with the prin- 
cipal material direction (1 or fiber direction) and X, Y and X 45 are the 
tensile strengths along, transverse to and at 45° to the principal 
material direction respectively. While equation (4) Is for uniaxial 
loading only, other stress states are considered by the author. 

Using assumptions similar to those of Ashkenazi, Hill developed the 
failure criterion for an orthotropic material under a generalized state 
of plane stress 
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where the stresses °]° 2 * T i? are along the material axes, S Is the In- 
plane shear strength and X and Y are as previously defined. 

A more general theory due to Puppo and Evensen can be expressed as 
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where y * -y- Is an Interaction factor and the other quantities are as 
previously defined. Further, the authors Indicated that their theory 
could be adapted to a wide range of materials by wrltlnq the Interaction 

f3s 2 r 


factor as y * . Equations (6) are for orthotropic materials but 

the general theory Is not so limited. 

The failure theories given above are used In connection with lami- 
nation theory in the following way. The properties of a single ply or 
lamina must first be determined. As single plies are quite thin, ^ 0.005 in, 
(0.127 mm), properties are usually determined by testlnq a multi-layered 
unidirectional laminate. These are taken as the assumed properties of a 
single ply and the properties, A^j, of general laminates are then calcu- 
lated for each configuration being considered. (Obviously, if the 
properties of the unidirectional laminate do not represent the properties 
of a single ply or If different defects are present in the manufacture 
of the general laminate as opposed to the unioi rectional laminate, the 
"calculated" properties are likely to be in error.) From a knowledge of 

O 

the experimental loads and the A^j, the in-plane strains, c j , can be 
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calculated from equation (3). These can be used to determine stress, 

, In the layer from equation (2) which can then be substituted Into 
the appropriate failure criteria. Thus, each failure theory Is then 
expressed In terms of the applied loads Nf and a general laminate can be 
examined ply-wise for each load Increment until a lamina reaches the 
failure state, (of and Nf above refer In each case to qlobal coordinates.) 

The dlffitul* question encountered after first ply failure (FPF)^ 

is how unloading occurs. Petite and Waddoups^ use the negative tangent 

modulus technique to unload the failed ply together with a maximum 

strain concept to uncouple lamina stiffnesses. An alternate approach 
12 

due to Sandhu Is similar except It uses a maximum strain energy concept 
to uncouple the failed lamina. Another simple technique cited by Jones 3 
is to uncouple the degraded lamina by assigning zero stiffnesses to the 
three off-diagonal terms of equations (2). The appropriate diagonal 
terms are zeroed using a maximum strain criteria when and If a diagonal 
strain Is exceeded. 

We have elected to use a different approach to account for the 
reduced laminate load carrying capacity duf to FPF. In any general 
laminate certain plies are likely to fail .rst. For example In a 
[0°/±30°/0°]5 tested in uniaxial tension in the 0°, FPF likely occurs in 
one of the ± 30° plies. Each such ply would not necessarily fail at the 
same time or have the same failure plane. It is assumed, however, that 
failure in each 30° ply would occur simultaneously but planes would be 
staggered as shown In Fig. la. Final failure is assumed to occur as 
shown in Fig. lb. The analysis procedure of this staggered FPF model 
follows. 
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Using lamination theory and ar.y of the previously mentioned failure 
criteria, the laminate Is loaded Incrementally until FPF occurs (e.g., 
the 30° plies for a [OV^OVO 0 ^ as shown 1r, Flq. la). The remote 
stress for uniaxial tension of a li nate prior to FPF Is given ^s 


P 

°n ' WT 


(7) 


where P, W and T are the applied load and the laminate width and thick- 
ness respectively. After FPF, the stiffnesses of the failed laminae 
(♦ 30°) are zeroed and the remote stress Is taken to be 


P_ 

WT' 


( 8 ) 


where T' Is the new assumed laminate thickness. For the [OV^OVO 0 ]^ 
example of Fig. la, T * 8t and T' * 6t where t Is the lamina thickness. 
Note that the new thickness T' ■ 6t assumes that the "failed" ± 30° plies 
are still capable of carrying a fraction of the applied stress and are 
not assumed to be completely Ineffective. The laminate continues to 
carry load until final fracture as shown In Fig. lb. 

The procedure described above can be used for any laminate. Of 
course, for a unidirectional laminate, FPF Is equivalent to total failure. 
Because Interlaminar stresses are neglected this approach gives slightly 
conservative results for qeneral laminates as will later be demon- 
strated. 

Nonlinear lamina material behavior is easily Incorporated in 
laminate failure predictions due to the Incremental nature of the compu- 
tational procedure. In addition, the procedure can be used to predict 
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the uniaxial stress-strain response of a laminate for the same reason. 

Experimental Procedures and Results 

The materials studied In this Investigation were manufactured by 
Lockheed (Sunnyvale, Calif.) from prepreg tapes composed of Hercules 
(Magna, Utah) graphite AS fibers and epoxy resin 3501. The fundamental 
properties of the fibers were: 380-400 ks 1 (^ 2620-2760 MPa) tensile 

strength, 30-40 x 10^ ksl (^ 207-276 MPa) elastic modulus and 10,000 fibers/ 
tow. The resin was a hot-melt 100?; solids epoxy. No properties of the 
resin were available. The laminae were about 65% fiber by volume. The 
resuiting [0°]gs and [OVOOVO 0 ^ laminates were medium strength-medium 
modulus composite materials, .arge 0.80 In (20.32 mm) thick plates were 
received from which Individual specimens were machined. 

Machined surfaces were examined by a scanning electron microscope 
(SEM) and various Inherent flaws were found. Figure 2 shows two such 
flaws which gives an Indication of the size of those found. 

Specimens of dimensions 8 In x 0.5 In x 0.08 In (203.2 mn x 12.7 ran 
x 2.032 ran) were cut from both the [0°]gs and [0°/* 30°/0°]2s panels with 
three different orientations. These were In the principal (0°) fiber 
direction, at 45° to the principal fiber direction and at 90° to the 
principal fiber direction. After machining, all specimens were stored 
In a desslccator until tested. Specimens were allowed to sit and stabilize 
to the test environment for at least 1 hr prior to testing. Test tempera- 
tures were generally at room temperature of approximately 75°F and the 
relative humidity was generally less than 60%. The specimens were 
Instrumented with longitudinal and transverse strain gages and were 
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tested In uniaxial tension at strain (head) rates from 0.002 In/mln (0.0508 
nri/mln) to 2 In/mm \b0.8 mm/mln). Both [0°]gs and [0°/ , 30°/0°] laminates 
were tested at the various rates and with loads at 0°, 45°, and 90° to the 
principal material direction. All tests were conducted without tabs, uslnq 
sandpaper between epoxy coated wedqe arlps to minimize penetration of grip 
serrations Into the graphite/epoxy materials. Data collected In this 
manner correlated well with data collected by NASA-Ames using tabs. 

Because of the variability of stress-strain data, a statistical 
subroutine from the library of our IBM-370 computer was used to condition 
the data. Axial loads and longitudinal and transverse strains were fed 
Into the computer program. The program calculated the axial stress, 
fitted the stress-strain data with a polynomial equation using a tech- 
nique similar to the least squares method, qave a listing of the poly- 
nomial coefficients and plotted the experimental data toqether with the 
fitted curve. Most data could be fitted well with a third order poly- 
nomial and only a few required a fourth order curve for a qood fit. The 
gradient of the various stress-strain (o x - e x ) or strain-strain 
(c x - ty) polynomials at any point was taken as the Instantaneous laminate 
modulus or Poisson's ratio respectively. These instantaneous values for 
the unidirectional laminate were used as input to lamination theory to 
obtain the response of the general laminates. 

Generally, the Instantaneous modulus, do/de * E t , and Poisson's 
ratio changed throughout each test by as much as 25%. While decreases 
usually occurred. Increases were noted for the [0°] 8S and [0°/* 30°/0°] 2s 


tests. 
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Final fracture surfaces were In general perpendicular to the load 
direction for the [0°]g S and [0°/ • 30 o /0 o ]«5 test as is shown In Figure 3. 
As may be seen from the same figure, longitudinal splitting was always 
at least a secondary (and sometimes a primary) fracture mode for the 
[°°]8S tests. All other fracture surfaces for other fiber orientations 
tended to follow the principal fiber direction. 

The fracture process for the [0°/ • 30°/0°]2 S specimens was often 
quite explosive and In many cases fragmentation occurred; l.e., fragments 
would be expelled considerable distances from the test machine. Also, 
for this series of tests, post specimen examination did reveal some evi- 
dence of delamination. However, It Is felt that this delamlnatlon 
occurred after separation due to rebound Impact forces (see Figure 3 
which she • specimen which fractured In two places simultaneously). 

In many tests various factors Indicated that fracture of Inner 
plies occurred first. Figure 4a Is an SEM photograph of the central por- 
tion of a [0°]g5 specimen fracture surface and flqure 4b Is an SEM photo- 
graph of Individual broken and pulled out fibers. 

For the [0°/±30°/0 o ] 2s , [45 °/ 1 5°/75°/45° D 2S , and [90°/ ♦ 6 n°/ 90 °] 2s 
tests audible noise could be heard In nearly every case long before separa 
tlon. Undoubtedly this audible noise was due to breakaqe of either Indi- 
vidual fibers or individual plies and most likely was due to the latter. 

Another indication that fracture occurred on Inner plies first can 
be seen by examination of Figure 5 which shows the load-time and strain- 
time traces from one strain-rate test for a [90 < 7‘60 0 /90 0 ],,^ specimen. As 
may be observed, small excursions in load occurred simultaneously with 
large excursions in strain. This was accompanied by large amounts of 
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audible noise with no visible signs of failure on the surface of the 
specimen. A possible explanation for these observat':. , Is that fracture 
may have occurred first on an Inner ply or plies near the location of 
the strain gage on a surface ply resulting In a transfer of strain from 
Inner to outer plies with little loss In load carrying capacity. Such a 
porcess could cause a large Increase In strain on an outer ply at the 
strain gage site and not affect the load carrying capacity appreciably. 

It should be noted that the manufacturer's specification for the 
tensile strenqth of the [0°]^ material was 240 ksl (v 1650 MPa). Our 
tests, substantiated by NASA-Ames, Indicated a tensile strength of only 
151 ksl ('v- 1050 MPa). This extreme difference seems to be a fault of 
the manufacturing technique. Also, these graphite/epoxy panels contained 
extreme thickness variations, l.e., as much as 35*,. Thickness variations 
of Individual specimens were considerably less, l.e., 'v- < 10^. Such 
variations in strength and thicknesses are obviously very Important in 
attempts to understand the failure and fracture behavior of composite 
materials. 

F1gur f "> 6 through 11 show the stress-strain curves obtained for 
the [0°] 8S , [45°] 8s , [90°] 8s , [0°/* 30°/0°] 2s , [45 0 /15°/75 0 /45°] 2S and 
[90°/t60°/90°] 2 5 axial tension specimens at the various head rates indi- 
cated. Each stress-strain curve was obtained by computer aided statistical 
conditioning of three separate sets of tension test data for each configu- 
ration tested as discussed briefly in the last section. Several interest- 
ing facts can be discerned from an examination of Figures 6-11. The 
curves for :he [0°] 8S and [n o /±30°/0 o ] 2 5 laminates of Figures 6 and 9 
indicate very little rate effect. The differences noted In these 
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figures were probably due to experimental scatter. When none of the 
fibers were In the load direction varying degrees of rate effects were 
apparently present. Figures 7 and 10 tend to Indicate that the [45°] 8S 

and C45 0 /15°/75 0 /45°32S materials First show an increase of properties 

with Increasing rate, 0.002 In/mln (0.0508 mm/mln) and 0.02 In/mln 
(0.508 mm/mln) followed by a decreaa. of properties with Increasing rate, 

0.2 In/mln (5.08 mm/mln) and 2.0 In/mln (50.8 rrr/nln). The [90 o /*60 o /90 o ] 2$ 
data of Figure 11 tend to show that properties essentially decrease with 
Increasing rate. Examination of Figure 8 for the [90°]g$ laminate Indi- 
cates. as expected, that when a laminate was matrix dominated properties 
In fact Increased with Increasing r«ite. Apparently some type of transition 
from a fiber dominated laminate with no rate effect to a matrix dominated 
laminate with rate effect did exist. The exact form of the transition 
Is unclear but It does appear that If more than 50? of the fibers were In 
the load direction little rate effect was present. Further, it would 

seem that for certain stacking sequences, i.e., [45 0 ]2$, [45 0 /75 0 /15 0 /45 0 32S 
and [90°/*60°/90°]2s that Increased rate had a deleterious effect on 
the material properties. 

Bilinearity Is evident in Figures 6 and 9 for the [0 °]qs and 
[0°/+30 3 /0°]2 s material where the stiffness tends to Increase with stress 
or strain level. This Is likely due to the srlssorlng effect of 
straightening the fibers or possibly the failure of the 30° plies in the 
[0°/+30 o /0°]2s case. Analytical predictions to be discussed tend to con- 
firm this latter observation. Inelastic stress-strain behavior was ob- 
served for the other orlentatons tested. Anotner interesting fact found 
was that Poisson's ratios were ntnative for the [45°/15 0 /75°/45 0 ]2s tests. ^ 
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Ana lysis Predictions Compared to T es t Resul ts 

As explained previously, an Incremental numerical approach was used 
to obtain the analytical stress-straln-fal lure predictions. The 
Instantaneous values of the elastic properties were obtained for a par- 
ticular stress-strain level from the polynomial representation of the 
stress-strain response of a unidirectional laminate. These values were 
substituted Into the classical lamination theory relationships (equations 
1 to 3) for the next Increment of the remote load. The resultinq stresses 
were substituted Into the respective failure criterion (equations 4 to 6), 
plywlse. When the state of stress In a ply reached or exceeded the 
Interaction surface predicted by the failure criterion, the ply was un- 
loaded either by the reduced stiffness method of Jones^ or the staogered 
fracture surface method described earlier. The laminates were assumed to 
have the same tensile and compressive properties. 

The results obtained by the above procedure are shown In Figures 
12-14. Figure 12 shows the comparison between analytical prediction* and 
experimental results for the [0°/*30°/0°]2s laminate. Excellent corre- 
lation Is seen up to about 50 ks 1 (344. C MPa). At this stress level, 
failure of the 4 30° plies was predicted and is Indicated by the steeper 
slope of the stress-strain curve. Thus, analysis and experiment indi- 
cated the same trend. After FPF, the analyses predicted higher stresses 
than were experimentally found. Reasonable ccrr.lation between theory and 
experiment was found for the [45 0 /15°/75°/45 o ]2s laminate as is seen in 
Figure 13. Excellent correlation for the [90°/i60°/90°]2s laminate was 
found as evidenced by Figure 14. On the whole, the analysis predicted 
larger stiffnesses than were experimentally found. This and the distinct 
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deviation between analysis and experiment after FPF for the [0°/‘30°/0°]2s 
laminate was ilkely due to omission of the Interlaminar shear stresses 
In the analysis. 

Both analytical unloadlnq techniques used gave reasonably the same 
stress-strain predictions. However, the staggered fracture surface method 
tended to give more conservative strenqth predictions than the reduced stiff- 
ness method as may be observed on examination of Table 1. This was primarily 
due to the simplicity of formulation of the former over the latter. The 
staggered fracture surface method Is, In general, computationally easier 
and a good bit more flexible when modeling laminate behavior beyond FPF. 

This would be true especially If interlaminar effects were Included. 

Among the three failure criteria used, Ashkenazi's appears to be 
the best. It should be noted, however, that equation (4) is only for 
uniaxial tension. The exponent of the interaction factor for equation 
(6) used herein was unity. By adjusting this exponent for the various 
laminates, identical analytical and experimental results could have been 
obtained. Unfortunately, for the laminates tested herein, this exponent 
was not a constant and thus the value of unity was used. However, the 
Puppo and Evensen theory appears to be most promising as it can be used 
for multl-axiil loading, allows for material variations through the 
interaction factor, and can be used on multi-axial laminates without a 
plywlse analysis. On the whole, good strength predictions were obtained 
from all three criteria using both unloading techniques. 

Conclusions 

SEM micrographs have been presented which indicate the types and 
sizes of inherent flaws encountered in the graphite/epoxy laminates 
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Investigated. Evidence has been presented which Indicates that first 
failure occurred on Inner plies. It Is thought that this Interior 
failure may account for the audible noise which could be heard at various 
stages of the loading process. As a result, perhaps accoustlc emission 
might represent an Ideal method for Identification of failure modes and 
mechanisms. 

Computer conditioned stress-strain data ha'e been presented which 
Indicated that fiber dominated laminated polymer matrix composites appear 
to have no rate effects whereas polymer matrix dominated laminated com- 
posites do show a rate dependence. Some type of transition from one 
behavior to the ether for other laminates does seem to occur. It appears 
that these laminates first show Increase of properties with Increased 
strain rate followed by a decrease of properties for additional Increases 
In strain rate. This Is similar to results found previously.^ Thus, 
high rates may have deleterious effects on some types of laminates. 

Three failure theories have been examined with two methods of un- 
loading. Classical lamination theory in conjunction with these theories 
has been shown to give reasonably good predictions to actual stress-strain 
response. Further, the theories also Indicated first ply failure on the 
Interior of the [ 0 °/+ 30 °/ 0 °] 2 s laminates with subsequent increased stiff- 
nesses. The three theories and the two unloading schemes have been shown 
to give reasonable failure strength predictions. It has been suggested 
that the Puppo-Evensen theory and the staggered failure surface method 
appear to be the better approach of those attempted. 
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Table 1. Predicted Failure Strength, Normalized with Respect to the 
Experimental Load. 


Laminate 

Orientation 

Ashkenazi 

Hill 

Puppo-Evensen 

Staggered Failure Surface Method 


[07^3070°]2S 

0.979 

0.900 

1.021 

[4571 5775745° ] 2 $ 

0.912 

0.879 

0.889 

[90°/±60790°] 2S 

0.636 

0.648 

0.587 


Reduced Stiffness Method 


[0°/±30 o /0°] 2s 

0.971 

0.965 

0.968 

[45°/15775°/45°] 2s 

1.002 

1.168 

1.054 


[907-*60790°] 2S 


0.827 


0.786 


0.864 
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Fig. 1. Sideview of [0°/* 30°/0‘]5 Laminate. 
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Figure 6. Stress-Strain Curves of [0°]g$ Laminate. 
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Figure 7. Stress-Strain Curves of [45°] ft r Laminate 
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Figure 8. Stress-Strain Curves of [90°]g5 Laminate. 
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Figure 9. Stress-Strain Curves of [0°/i30°/0°]2s Laminate. 
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Figure 10. Stress-Strain Curves of [45 0 /15 0 /75 0 /45 0 ]2s Laminate. 
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Figure 12. Comparison of Experimental Results with Analytical Predictions 
of [0 o /±30 o / 0°]25 Laminate. 
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Figure 13. Comparison of Experimental Results with Analytical Predictions 
of [45 0 /15 0 /75°/45 0 ] 2s Laminate. 
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Figure 14. Comparison of Experimental Results with Analytical Predictions 
of [90 o /*60 790 o J 2s Laminate. 
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